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Abstract 

The average density of states (DoS) of the one-dimensional Dirac Hamiltonian with a 
random mass on a finite interval [0, L] is derived. Our method relies on the eigenvalues dis- 
tributions (extreme value statistics problem) which are obtained explicitly. The well-known 
Dyson singularity g{e; L) ^ ~ j^|i^T]7y is recovered above the crossover energy ~ exp — 

Below Ec we find a log-normal suppression of the average DoS Q{e; L) ~ j7j^/x ^■'^P(~r ^^^'^ I^D' 
PACS numbers : 72.15.Rn ; 73.20.Fz ; 02.50.-r. 

1 Introduction 

Let us consider the one-dimensional Dirac equation WdV' = ^4' fo^' the Hamiltonian 

Hb = -a^idx + P(t>{x) , (1) 

acting on a two-component spinor. (f){x) plays the role of a mass and will be taken as ran- 
dom. The Dirac matrices are chosen as /3 = 7*^ = ui and ax = 'y^J^ = —(^2, where cjj are 
the usual Pauli matrices, corresponding to the representation of the Clifford algebra 7*^ = ai 
and 7^ = — i(T3. One-dimensional random Dirac Hamiltonians appear in several contexts of con- 
densed matter physics ranging from disordered half-filled metals [H |2] , random spin-chain mod- 
els (random antiferromagnetic spin- 1/2 chains |3lHl|5], random transverse field Ising spin- 1/2 
chains jH [3 18] , spin-Peierls chains and spin-ladders [HI [10] ) , or organic conductors lHj . Another 
common representation of Dirac matrices /? = di and ax = crs, related to the one chosen in the 
present article by a unitary transformation lA = {1 + i(Ti)/\/2, shows that the Dirac equation 
has the form of the linearised Bogoliubov-de Gennes equation describing a superconductor with 
random gap j l2| [T3] . We obtain a connection to another important, and well-studied problem 
by squaring the Dirac Hamiltonian: "H^ = —9^ -|- (/)(x)^ -|- \-)^(f)'{x) = —d^ + (j){x)'^ + (T30'(x). 
This leads to a couple of isospectral supersymmetric Hamiltonians H± = —-^ + 4>{x)'^ it 0'(x). 
The corresponding Schrodinger equation can be transformed into a Fokker-Planck equation 
dtP = dx{dx =F 2(j){x))P describing classical diffusion in random force field, studied in numerous 



works |14| [TE\ [THl [T71 [TE\ El dSl EO] . If the random mass is a Gaussian white noise with charac- 
teristics (j){x) = fig and (l){x)(j){x') — (j)"^ = g 5{x — x'), where denotes the average with respect 
to its reahsations, an exact analytical expression for the density of states (DoS) was obtained 
da HZ] (see §8.2 of Ref. [13,) [] Another solvable case is the one where the mass is chosen 
as a telegraph noise [22] [i.e. with exponentially decaying correlations). The same low energy 
properties are obtained. This is related to the fact that any short-range correlated noise, upon 
large scale renormalization, reduces to a Gaussian white noise. 

Here we focus on the case of Gaussian white noise with (t){x) = which is known to yield a 
Dyson singularity for the DoS Q{e — > 0) ~ |e| in^g/kl) ^^"^ ^ delocalisation transition at e ^ 0. 
The aim of the present article is to study how the Dyson singularity of the DoS is affected by 
boundary conditions. We denote by f?(e; L) the DoS of the Hamiltonian ([T]) on a finite interval 
[0, L\. We will obtain the behaviour of the average DoS Q{e] L) for a finite length. 



2 Boundary conditions 

Let us specify the two different types of boundary conditions we shall consider. 

Type (D) .-A simple way to introduce boundary conditions ensuring confinement in a domain 
T) is to consider the Dirac equation [i^ — (j){x)]tp{x,t) = with infinite mass 4'{x) = <I> — > oo for 
X outside T). This leads to the so-called "bag model" of hadronic physics [23 : the boundary 
conditions are {1 + m ■ ^)iP\qx> = 0, where ft is the unit vector normal at the boundary dT> 
of the domain. They force the vanishing of the component of the current density = ip^ip 
perpendicular to the boundary, where = Tp^j^. In our case, for x G M~ the stationary solution 
of the Dirac equation for a constant mass is a plane wave •il){x) = ue^^^ where the spinor u is 
solution of {oxk + f3^)u = eu, for e'^ = k"^ + In the limit $ ^ oo it becomes an evanescent 
wave with k ~ — i<I> ; this yields (1 — i'y^)ip{0) = 0. Similarly, at the other boundary we find 
(1 + i'y^)ip{L) = 0. We denote these constraints as "type (D) boundary conditions". With our 
representation of the Clifford algebra, they read (1 — a3)^p{0) = (1 -|- a^)'4j[L) = 0, i.e. each 
of the two components of the bispinor tp = (ip, x) must vanish at one side : x(0) = '^{^) — 0- 
This condition obviously ensures the absence of Dirac current flow = —i/j^a2'ip across the 
boundaries. 

Type (S) : A. second interesting choice is given by (1-|-(T3)'0(O) = {l + a^)'4j{L) = 0. These condi- 
tions coincide with Dirichlet boundary conditions for the associated supersymmetric Schrodinger 
Hamiltonian H± and will be denoted as "type (S)"[^ 

^ Note also the Ref. [2T] where the most general Id Dirac Hamiltonian Hd = 02 i^ + ai4>{x) + cr3A(a;) + V{x) 
was studied for <j), X and V Gaussian white noises. This study showed that there is no other DoS singularity than 
the one arising in the random mass case. 

^ General boundary conditions : It is possible to set up more general conditions forcing the absence of proba- 
bility current flow at the boundaries : these are the set of conditions parametrised by some real number i9 [24] : 
(e"''' + in ■ 7)i/'|s-D = where 7^ is the chirality matrix, given by 7^ = 7''7"'^ in dimension d = 1 + 1. These 
general boundary conditions are obtained as follows : one considers the Dirac equation in a bounded domain 
T). The boundary conditions can be written in • "ftpld-n = {A -\- B y^)-tp\gT>, where n is the unit vector normal 
to the domain, A and B two complex numbers. The fact that (n • 7)'^ = —1 leads to A'^ — B^ = 1. Upon 
imposing the normal component of the Dirac current to vanish n • J7^|ax> — 0, we find A G M and B G iM, and 
thus the desired form. With the representation chosen in our article these general boundary conditions take the 
form (e-'''±'^2 ± a3)tpix±) = with X- = and x+ = L. For type (D) we have '&+='&-= 0, and for (S) 

— TT, i9+ — 0. Since j?± 7^ or tt breaks the particle-hole symmetry, this case will not be considered here. 
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3 Eigenvalue distributions 



We now consider the Dirac Hamiltonian ([T]) on [0, L] for the two kinds of boundary conditions 
introduced above. The particle- hole symmetry takes the form ip ^ ip = a^ip : a^TLna^ = — T^d- 
Since both boundary conditions preserve the particle-hole symmetry, eigenvalues appear in pairs 
(by convention we choose n = 1, 2, 3, • • • and > 0). We denote by ^(e) the DoS per 
unit length for an infinite volume. It is related to the DoS of the finite-size system by £»(e) = 
lim^^oo L)/ L. Because of particle-hole symmetry we can restrict ourselves to e > 0. 

Our method to obtain the average DoS Q{e] L) relies on the evaluation of the eigenvalue 
distributions W„(e;L) = (5(e — €^[0, L]), derived in Ref. [25] for (S)-boundaries. Finding the 
distributions of the (ordered) variables corresponds to an "extreme value problem" (here 
for correlated variables). We recall the idea of the method : let us imagine that we impose 
the boundary conditions to hold solely at one end of the interval, say x = 0. As x increases 
the two components of the spinor ip{x) = (y?(x), x(^))) which solves the Dirac equation, vanish 
alternately ; we denote by the distance between a node of ip{x) and the closest node of x(^)- 
Since the evolution of ip{x) is Markovian, the A^'s are identical and independently distributed 
{i.i.d.) random variables whose statistical properties are obtained by solving a first passage time 
problem (see Ref. [55] for details). Let us introduce the distribution ttm^v) of the sum of M 
i.i.d. rescaled lengths y = (Ai -|- • • • -|- Ajv/)AA(e), where 

is the integrated DoS per unit length for an infinite volume ; C = 0.577... is the Euler-Mascheroni 
constant (the exact expression for AA(e) may be found in Refs. |12| [TS] [T7]). This rescaling is 
motivated by the fact that Am = l/(2AA(e)) [25j, and hence y = M/2. Both boundary conditions 
at X = 0, L are satisfied whenever the sum Ai -|- • • • -|- Am coincides with the length L. Therefore 
the distribution yVn(e; L) is given by|^ 



Wnie;L)=LQ{e) Wn[LMie)j , (3) 

where the dimensionless function Wn{y) is related to the distributions TTMiu) as follows : 

Type (S) : The same component of the Dirac spinor must vanish at the two sides of the interval. 
Therefore L must coincide with a sum of an even number of lengths A^'s (see eq. (132) of 
Ref. t25j) : 

tn(^)(y)=vr2n(2/)|. (4) 



Type (D) : Both spinor components must vanish at one side of the interval. Hence L must 
coincide with a sum of an odd number of lengths A^'s : 



n7f)(y)=vr2„_i(2/) . (5) 



The distributions '7r„(y) are explicitly known in the low energy limit e <^ g, i.e. for gL ^ 1 
where Wn(e; L) is essentially concentrated below g. Indeed, the characteristic function of the 
lengths Am's is given by [2^ e~"^ ~ cosh""^ y^a/J\f{e). Whence we may write 

f dq e'^y 
" 2i^ 

^ Note that the change of variable from e to LM{€) corresponds to the so-caUed spectrum unfolding leading 
to a unit density of variables. 
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where the integration is taken along the Bromwich contour. Starting from this integral repre- 
sentation we obtain the following explicit formulae 



-,2k 



n + k \ ("+1/2)^ 



(7a) 




(7b) 



whose demonstration is provided in appendix [A} 



4 DoS on a finite interval 

The average DoS can be related to the distributions recalled above : 



Q{e;L) = f2^nie;L). 



n=l 



Summation over n leads to 



g{e;L) = Lg{e) v[LM{e)) 



(8) 



(9) 



where the dimensionless function T>{y) depends on boundary conditions. Note that this DoS 
has some interest for studying the problem of Id classical diffusion in a random force field 
with dilute absorbers [261 (see also Ref. [27] for an analysis of this problem with the real space 
renormalisation group) . 

Boundary conditions of type (S).- The summation of the (S)-type distributions (|4| gives 



B 2m sinh ^/g 



(10) 



The integrand is meromorphic in the complex plane (there is no branch cut since yfq is the 
argument of an even function). The integral can be computed from the residue theorem. The 
integrand possesses a single pole at g = with residue equal to unity and an infinite number of 
double poles on the real axis at q = = — (nvr)^, with n = 1, 2, 3, ■ 



^{q — qn)"^, we find the residues Res [ 



4(3„ 

fore we have 



sinh 



dq I sinh^ ^ J q=gn 



Using that sinh ^ ~ 

q^qn 



Vs{y) = 1 + 2 ^ (1 - 2(n^)2y) e-("-)'2' 



Try 



= (2+4g„)e''"S'. There- 

3/2 E (11) 



n=l ^ ^ n=l 

where the second series expansion may readily be found from Poisson's summation formula (see 
appendix |B]) . We may also check that the summation of ( 7b ) leads to (11). 

Boundary conditions of type (D).- Compared to the (S) case, an additional coshy^ appears in 
the integrand upon summation of ^ : 



dq cosh ^ 

— 2^ I 

S 2m sinh y/q 



(12) 
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Figure 1: Functions T>s{y) (black continuous line) and 'D£,{y) (dashed blue line). First distri- 
butions Wn{y) are also plotted in thin lines. 



It adds a sign (—1)'^ to the residues obtained in the previous case and we thus find 



Voiy) = 1 + 2 ^(-1)" (1 - 2(n^)2y) e-("-)'2' 



n=l 



^ n=0 



(13) 



The two different boundary conditions can be treated on the same footing by writing 



2?5,D(y) = (l + 2y-^) j;(±l 



(14) 



and 



Try 



3/2 



with r] 



n=0 



for S 
1/2 for L» 



(15) 



This second series expansion allows to analyse the low energy behaviour. It is worth noticing 
that this representation for 2?(y) is very similar to the corresponding one for vji^y), for both 



kind of boundary conditions. In the (S) case eqs. (22) and r.h.s. of (13) only differ by the 
(-1)"+^ in the sum, while in the (D) case ^^ and r.h.s. of ([20]) differ by a (-l)"(2n + 1). As 



a result we have 'D{y — > 0) ~ '^i{y) in both cases. The two functions are plotted on figure [1} 
'Ds{y) presents a monotonous behaviour what leads to a diminution of the low energy DoS for 
small energies. Interestingly for (D) boundary conditions, while exponentially suppressed for 
y — 0, 'D£){y) increases for intermediate values of y ; as a result (D) boundary conditions induce 
an increase of the DoS at intermediate energies. 

The low energy DoS presents a log-normal suppression : 



and 





16 p-^ 1-^(9/1^1) 
e \j2TTgL 




|e| = 5 6-^9^2 


Ssie; L) ~ 


for 






e y/2-KgL 




|e| ^e'^ = ge-^ 


to(e; L) ~ 


for 



(16) 



(17) 



The DoS reaches its maximal value at e* the (S) case and at ~ ge in the (D) 

case. 

The weaker log- normal suppression in the (D) case as compared to the (S) case is due to 
the fact that for a given realisation of the disorder, the spectra for the two different kinds of 
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boundary conditions are such that ef < ef < < ef < • • • , where {en'^} denotes the spectrum 
for boundary conditions of type (D) and (S) respectively. 

Since T){y oo) ~ 1, eqs. ( 11|13 ), we recover the Dyson singularity for intermediate energies, 
as expected : 

Qie; L) ~ L Q{e) ^ for « |e| « g . (18) 



At higher energies |e| ^ one should recover the free DoS : g^e; L) ~ L/vr. 

The average DoS £»(e; L) is represented for various values of L on figure [ij For the lowest 
energies e < e* ~ ge~^^ the DoS is suppressed in both cases. In the intermediate range e* < e < 
ec ~ ge~^^, the effect of the boundaries is to reduce the DoS in the (S) case but surprisingly 
to increase the DoS in the (D) case. 

As L is decreased the Dyson singularity is rapidly converted to a strong depletion of the low 
energy DoS. For (S) boundary conditions this occurs for a surprisingly relatively large length 

~ 10/ g. For the (D) case the memory of the Dyson singularity persists up to smaller lengths 
since an increase of the low energy DoS is apparent up to L'^ ^ 1/g (figure |2]). 
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Figure 2: Average DoS per unit length for a finite L ; parameter is g = 1. Left : (S) boundary 
conditions for sizes L = oo, 30, 25, 20, 15, 10 & 5. Right : (D) boundary conditions for sizes 
L = oo, 12, 6, 5, 4, 3, 2 & 1. 



5 Conclusion 

By using known results on the energy-levels distributions W„(e; L) we have obtained the average 
density of states of a random Dirac Hamiltonian on a bounded domain [0,L]. For gL ^ 1, we 
have seen that the DoS does not feel the boundary as energies become larger than the energy 
scale ec ~ ge~^^. In the intermediate range < e < ec, where e* ~ ge~^^, the DoS may 
be diminished by boundaries, (S) case, or increased, (D) case. For the lowest energies the DoS 
presents a log-normal suppression for both kind of boundary conditions. 

Note that the effect of one boundary condition was already studied in Ref. [10] where the 
average local DoS on a semi-infinite line was computed by the Berezinskii technique. These 
authors found an increase of the local DoS by a factor ln(5f/|e|) S> 1 close to the boundary (at 
a distance much smaller than the scale g~^ln'^{g/\e\)). Here, by imposing boundary conditions 
at the two sides of finite interval, we have shown that boundary conditions can induce both an 
increase or a decrease of the DoS. 

It is interesting to interpret our result within the classical diffusion problem. The (S) bound- 
ary corresponds to diffusion in a random force field with two absorbing boundaries while the 
(D) case corresponds to one reflecting boundary and one absorbing. The return probability 
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is the Laplace transform of the DoS, therefore the log-normal suppression of the average DoS 

coincides to a log-normal decay of the average return probability Jq dx P{x,t ^ oo\x,0) ~ 
exp-^ln2(52t). 

To close this paper, let us examine the relation between our main result 
localisation properties. It is only below the crossover energy Cc ~ ge~^^^ that eigenstates do 
feel the boundaries. The transition corresponds to the case where the localisation length is of 
order of the system size ~ -^i therefore this criterion suggests that the eigenstate of energy e 
is localised on a scale ~ g~^ln^{g/\e\). 

The question of localisation was studied in several works (Refs. 128\ [29] for a tight bind- 
ing Hamiltonian with random hoppings and in Ref. ^/7j for the continuous supersymmetric 
Hamiltonian) that have demonstrated the vanishing of the Lyapunov exponent |^ at the band 
edge 7(e ^ 0) ~ 5 ln~^((5r/|e|). The Lyapunov exponent provides another possible definition 
of the inverse localisation length == I/7, therefore much shorter than the previous one : 
~ g^^ ln(g/\e\) <^ In other terms the Lyapunov exponent analysis would put the "locali- 
sation threshold"!^ at e* ~ ge~^^, much below than ec ~ ge~^^^. 

The existence of two characteristic length scales was pointed out in several works [T71 [3l [5l [10] . 
These authors have found that the average Green's function (i.e. the two point correlation 
function) decays over|^^e ~ ln^((7/|e|). The existence of two length scales was interpreted 
by Fisher I^SJ as a consequence of fluctuations (see also [5]). The "typical" localisation length 
characterises the decay of a typical wave function (more precisely is related to the average of 
the logarithm of the wave function) and the "average" localisation length controls the decay 
of the average correlation functions. Since we have considered an average quantity g{e; L), the 
fact that we have extracted the scale is consistent with Fisher's argument. 

Several arguments support the existence of a delocalisation transition at e — > 0. Yet they 
are of quite different nature and it is interesting to provide a brief review. The first four 
arguments are bulk properties of the model, (i) The Lyapunov exponent vanishes at low energies: 
7(e) = ^ for e — > [T7]. (ii) The calculation of the average Green's function [T71 [5] [TO] . 
{in) The DC conductivity of the model was computed at the Dirac point in Ref. ^ [2] and was 
found to be finit^ (iv) The statistical properties of the zero mode wave function [331 El IMj 
indicate long range power law correlations. Another two arguments are obtained from scattering 
analysis : {v) The distribution of the transmission probability through a finite slab of length 
L at zero energy was obtained in Ref. ^5] ; in particular the average transmission decreases 
like that is slower than the behaviour 1/L for a quasi-ld weakly-disordered conducting 

wire, {vi) The time delay distribution presents a log-normal distribution at zero energy |35| [36] 
Pl{t — > 00) ~ -^^^^ ^^P ~8gL ^^'^(9^)^ with moments diverging with the length of the disordered 
region |8| [36] t(0) = 2L and r(0)" ~ g "e^*^ ^ . Finally, two arguments consider the problem 

*We recall that the Lyapunov exponent characterises the exponential growth rate of the envelope of the wave 
function |3UI I13j 7 =^ Mmx^oo where H controls the envelope of the wave function tpix) — 6"'^' x oscillations. 

Note that in the scaling theory of localisation [31] the notion of a localisation threshold (or mobility edge) 
designates the energy separating localised and delocalised states for the infinite system. Within this usual termi- 
nology the model we are studying is always in a localised phase apart strictly at e = 0. What we denote here by 
"localisation threshold" is the energy separating localised and delocalised states for a finite system size L. With 
this definition there can not be a sharp localisation transition, due to fluctuations over disorder realisations. 

^Note that a related quantity first obtained Refs. [15II16| is the propagator, the inverse Laplace transform of the 
Green's function, or more precisely the average conditional probability P{x,t\x' ,0) = (j)o{x)<l)o{x')^'^{x |e^*^+ 1 x' ) 
of the related Fokker-Planck equation, where (j}o{x) is the zero mode of the supersymmetric Hamiltonian H+ ; the 
result of these works was reproduced in Ref. |32] from stochastic Riccati analysis and in Ref. [19] with the real 
space renormalisation group. 

In the localised phase in Id, the AC conductivity vanishes at small frequency according to Mott's law 
ReCT(cj) ~ tj^ln^w [30] , 

® This suggests that a particle of energy e « injected in a disordered region of size L will spread over the full 
interval and remain trapped a very long time depending on L ; whereas in a strongly localised phase the particle 



9fl5) with the 
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on a finite interval taking into account boundary conditions : {vii) The study of extreme value 
statistics of energy levels reveals spectral correlations for e — > [25] , whereas a localised phase 
is characterised by absence of level correlations [38]. {viii) Finally the present work has shown 
the influence of the boundary on g^e 0; L) due to delocalisation. 

It remains an interesting issue to go beyond this analysis of the localisation characterised 
by the two lengths and and comprehend better the role of fluctuations, for example for 
the DoS g{e;L). Another challenging issue would be to extend our results to two dimensions. 
Since our method is specific to one-dimensional systems, such an extansion would necessitate 
the development of other methods. This might be of interest due to the recent revival of two- 
dimensional random Dirac Hamiltonian physics, motivated by the study of Graphene. 
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A Extreme value spectral statistics 

We give the first distributions introduced in the text. We compute the integral ^ using the 
residue theorem. This requires the expansion of the denominator near the pole k„ = — ^(2n -|- 

l)^ n e N : cosh = - Kn) [1 + S + ^ (^) V i^^^ 

algebra gives : 

°° 2 



n=0 



7ri(y) = ^(-l)"7r(2n + 1) e-'^^^n+i) y (^g) 
) 

y (-l)"(2n + 1) e-(2«+i)V4. (20) 

vr2(y) = y [7r\2n + ify - 2] (2n+i)^y (21) 

n=0 

A °° 

= ^E(-l)"'''"'^""'^' (22) 



^ n=0 



n=0 



2 y'-'y+-2 



i{2n+l)^y (23) 



vr3(y) = E(-l)"^(2n + l) 

^ y (-1)"+H2n + 1) [(2n + 1)^ - l] e'^^^+D^^ (24) 



vr4 



2^y3/2 

^ ^ n=0 

, , ^r^^(2n+l)% 2. x2 2 /7r2(2n+l)2 \ 4 

(y) = E — -y - 2^ + ^) y + 2 ( ^—o — - + 1 ) 2/ - 3 

n=0 L V / 



e-^(2n+i)2, (25) 



^ ^ n=0 

would only explore a typical size what is reflected in the fact that the time delay distribution reaches a limit 
distribution limL-^oa Pl{t) — -Poo(t) WI]. 
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(the distribution vri (y) was already obtained in Ref . [T] where the energy ei is interpreted as the 
gap of a spin chain. It is also related to the distribution of the random bond lengths in the real 
space renormalisation group procedure [21 El IISl EZ] ; ^2(y) and TTi{y) are given in Ref. 

Expressions ( |20|22|24|26] ) are obtained from eqs. ( |19|21|23|25| ) by using the Poisson for- 
mulae given in the next appendix and iT3{y) = {-2y^-^ - 3y + 2)'^i{y), T^2{y) = -2(2y^ + 

1) E:Loe-"^(^"+^)'^ and .,{y) = {ly% + S{y^ - |)|^ + 2, - |) T.n^o^-'^^'''^'^'' ■ 



A.l Generating function 

We propose now a method that allows for a systematic determination of the distributions vr„(y). 
Let us introduce the generating function 



Q{z,y) 



oo 

E 

n=l 



^qy 



jg 2m cosh y/q 



(27) 



where summation was performed in the convergence radius \z\ < 1. The function arccos is 
single valued in the convergence disk, therefore we can write the poles of the integrand as 
Qn = — (nvr + arccos z)^ with n € N. We compute the residues by using ]^ sinh ^/q^ = 

(— 1)";77 — ^-^"^^ — r (It Is helpful to notice that arccosfx it iO^) = arccos(x) =F 10"*" for x real in 

^ ' 2(n7r+arccos \ f v / v / i 

[—1,-1-1]), whence 



(nTT+arccos z)^y 



y dz 



(28) 



n=0 



In order to apply the Poisson summation formula (32) and get the generating function of 
the distributions with odd indices, we exploit the symmetry with respect to a change in sign of 
the argument ^^^^^(-i)" e-('"'^+^'^^^°^^)'2^ = -Y^-l_^^(^-i)n Q-{n-'^+^rccos(-z)fy . 



0{z,y) 



G{z,y) - Q{-z,y) 



z d 
2y dz 



-co v 

+00 



(riTT+arccos z)^y 



z d 
y dz 



/Try 



(n + l/2) 

(z) e y 



oo 
2 



n=0 



Here the r„(z) are the Chebychev polynomials of the first kind T„(x) 
may be rewritten as 

[n/2\ 



Tn{x) 



n 
2 



EM) 



fe=0 



k\{n-2k)\ 



(2x 



(29) 

cos(n arccos(a;)). They 
(30) 



where \x\ is the integer part. Upon insertion into (29) we obtain 



;0(^,2/) = E 



Vr2n+l(y) 



Using that 



n=0 



n 



Try 



3/2 E 

n=0 



(_l)«(2n + l)e" 



(n + l/2)^ 



E(-i)* 

fe=0 



n + k 
n — k 



(2z) 



2k 



(31) 



P 



for p > n, we may relax the constraint on the summation with respect 



to k and extract (7a). Note that eq. (7a) allows to recover (20) and (24). 

In order to evaluate vrn(y) with n = 2k recall that we deal with a distribution of n positive 
i.i.d. random variables. It follows that we may obtain vr2fc(y) by convolution of T^2k-i{y) and 
^i(y) '■ '^2k{y) = /(f dx 7r2/fc_i(x)7ri(y — x). The integration is fairly cumbersome but one may 
verify that it yields (7b). This completes the computation of iTniy) for all positive integers n 
as announced in the main text. One may verify that the summation of all Tr2k{y) (or 'K2k+i{y)) 
yields the densities given by eq. (15). 
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B Two useful Poisson formulae 



Let us start by recalling the well-known Poisson formula J2neZ fi''^) ~ X^neZ /l^^'^) for any 
function /(x) defined on M, with jf{k) = f^dx e~^'^^ f{x) its Fourier transform. Applying this 
formula we obtain 

1 ( — 

used for iTn{y) with even indices (set r] = and a = 1/2) and for the DoS (set rj = oi ij = 1/2 
and a = 0). For the distributions vr„(y) with odd indices we need (with tj = a = 1/2) 

y (n + a) e^i™'' e-"'("+°)'2' = — — ^ V (n - ry) e2i-(«-'?)" e"^^ . (33) 
'^-^ i(7ry)3/^ 



C A probabilistic interpretation of the result for type (S) bound- 
ary conditions 

It is worth pointing that the function T>s{y) can be interpreted as the integrated distribution of 
the maximum of a Brownian excursion. Let us denote by {x(t), ^ t ^ 1) such an excursion. 
We establish a relation between the distribution of the maximum of such an excursion and the 
function Vs{y). 



C.l Maximal height of a Brownian excursion 

Let us consider the distribution of the maximum M of a Brownian excursion (x(t), ^ t ^ 1) 
(a Brownian bridge constraint to be positive). It can be written as a ratio of two path integrals : 

""^^^ Vx{t) e-5/0 dti' JJ^(^(t))^(M - x{t)) 

J x(0)=xn j-_n 

Proba [x{t) ^ M] = lim —— — , (34) 

xo^O+ rx{l)=xo 1 1 

Jx(0)=xo 

where xq > is a regulator. The ratio of path integrals may be rewritten as 

Proba \x{t) ^ M] = lim ^^oj^ '^'l^o) (35) 

where Hq = acts on functions defined on IR+ satisfying Dirichlet boundary condition at 

X = and Hi = acts on functions defined on [0, M] with Dirichlet boundary conditions. 

It follows that IM 



Proba [x(t) ^ M] = (^^) ^n^e-Hif) . (36) 



n=l 



C.2 First-exit time 

We consider a Brownian motion on [0, l/\/2] with a reflecting boundary condition at x = l/\/2- 
We denote by r the time at which x{t) hits x = for the first time, starting from x = l/-v/2, 
and Tx the time needed to reach x = 0, starting from x. /i(x, q) = e~'J'^== obeys the BFPE 
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^■^h{x,q) = qh{x,q) with boundary conditions h{0,q) = 1 and dxh{x,q)\^^^ = [4Tj. We 
easily find : 

cosh - 1/^/2) 
^("''^) = ^^q ' ^^'^ 

therefore distribution of time r is given by inverse Laplace transform of 

1 



cosh y/q 



(38) 



We introduce the sum of n i.i.d such variables : y = ti + - ■ ■+Tn- The distribution of this variable 
was introduced in the text 7r„(y) = J^^^ ^ coS"^ ' where we have shown that 'Ds{y) = 

Yl^=i'^'in{y) = 0p^3/2 Yl^=i "-^ e~"^/^. Comparing with (36) we seen that the sum of these 
distributions coincides with the cumulative distribution of the Brownian excursion |£] 



Vsiy) = Proba x{t) M = TTy/yfl . 
It would be interesting to know whether this remark is purely accidental or not. 



(39) 
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